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Abstract 

We report on a theoretical study of quantum charge transport in atomistic models of silicon 
nanowires with surface roughness-based disorder. Depending on the nanowires features (length, 
roughness profile) various conduction regimes are explored numerically by using efficient real space 
order N computational approaches of both Kubo-Greenwood and Landauer-Btittiker transport 
frameworks. Quantitative estimations of the elastic mean free paths, charge mobilities and lo- 
calization lengths are performed as a function of the correlation length of the surface roughness 
disorder. The obtained values for charge mobilities well compare with the experimental estimates 
of the most performant undoped nanowires. Further the limitations of the Thouless relationship 
between the mean free path and the localization length are outlined. 
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INTRODUCTION 



Semiconducting silicon nanowires (SiNWs) are currently the subject of intense studies 
due to their prominent role in the downscaling limits of MOSFET devices, and also because 
they provide alternative materials to challenge quantum effects in low dimensionality 
Compared to classical planar technology, nanowires can better accommodate "all-around" 
gates improving field effect efficiency and device performances |2J. Vapor-liquid-solid (VLS) 
growth techniques have recently produced SiNWs with well controlled structural features 
[3], which opens up innovative approaches to the design of silicon-based nanodevices {4]. 



n 



Demonstrations of p — n junction diodes |5|, logic gates |6|, field effect transistors and 
nanosensors 8| have been reported. 

However, one key issue in the engineering of performant SiNW-based field effect tran- 
sistors (SiNW-FETs) is to ascertain how sensitive the charge mobilities are to structural 
features such as diameter, growth direction, and disorder. Surface roughness disorder (SRD) 



is a well known limiting factor in lithographic SiNW-FETs 



and its impact on ballis- 



tic transport in VLB-grown nanowires is a challenging awl important question Besides. 



0, 



12|. 



SRD effects also raise fundamental questions in the framework of localization theory 

Recent ab initio studies Q, [l4] have reported on specific surface effects, such as dopant 
segregation in small diameter SiNWs. However, these studies generally hardly cope with the 
analysis of the fundamental transport length scales in long disordered nanowires. Several 



theoretical works have also investigated the role of effective and simplified sur 
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ace disorder 



16|. 



models on the transport properties of nanowires-based materials or devices 

In this work, we report on a quantitative analysis of the transport length scales in atom- 
istic models of rough SiNWs. The description of the SiNWs is based on an accurate tight- 
binding Hamiltonian, previously validated by ab initio calculations 17| • The quantum trans- 
port properties are calculated with two different approaches. First, the elastic mean free 
path and the charge mobility are computed with an optimized, real-space order N Kubo- 
Greenwood approach [l8,[l^]. Additionally, the scaling properties of the Landauer-Buttiker 
conductance are investigated with a standard recursive Green's function method, to assess 
the effects of quantum interferences driving to the localization regime. Both approaches give 
complementary results and allow to explore a broad range of conduction mechanisms, from 
the ballistic to the diffusive and strongly localized regimes. 
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METHOD 



Description of the Surface Roughness Profile 



The SiNW Hamiltonian is a third nearest neighbor three center orthogonal sp 3 tight- 
rinding model that well describes the electronic structure of ideal (disorder-free) nanowires 



17l | . The SRD profile is defined as a random fluctuation of the radius of the nanowire around 



its average value Rq, characterized by a Lorentzian auto-correlation function, with a single 
intrinsic length scale [lo[ [^(J . In cylindrical coordinates, 



where: 
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^Pnk £ [0, 27r[ is a random number, L is the length of the nanowire and L r is the correlation 
length of the SRD. The silicon atoms outside the envelope defined by Eq. (jla|) are excluded 
from the nanowire and the dangling bonds are saturated with hydrogen atoms [2l[. In the 
following, we set Rq = 1 nm and renormalize the a n fc's so that (5R 2 (z,8)) 1 ^ 2 = 1 A, leaving 
L r as the only free parameter. The effects of the SRD on the transport properties have been 
analyzed for [110] oriented SiNWs with L r 's ranging from 0.54 to 4.34 nm [Fig. Ufa)]. 



Transport methodologies 

In the semiclassical transport theory, disorder effects can be characterized by the scatter- 
ing rate between the eigenstates of the ideal system [221] . The scattering time r is usually 
computed with the Fermi's golden rule, and by virtue of Matthiessen's law, can be split into 
an elastic (r e ) plus an inelastic (rj) contribution (1/r = l/r e + 1/Tj). The SRD is expected 
to dominate backscattering at low temperatures jsj], while the inelastic electron-phonon 
coupling plays a major role at room temperature. In the Kubo-Greenwood approach, the 
scattering time r e (E) and the mean free path £ e (E) = v(E)r e (E) are extracted from the 
saturation of the quantum diffusivity D(E,t) = AZ 2 (E,t)/t — ► 2v 2 (E)r e (E), where v(E) is 
the average velocity and AZ 2 (E,t) the quadratic spreading of wave packets with energy E 



3 




FIG. 1: (color online), (a) SRD profile and isodensity surface of the ground-state conduction band 
wave function of small sections of [110] SiNWs, for two values of L r . (b) Corresponding conduction 
band DoS (dashed and solid lines). The DoS of the ideal nanowire with radius Rq = 1 nm is also 



shown (dotted line), (c) Landauer-Biittiker conductance of the same ideal nanowire 
with the number of conducting channels (staircase line) . 
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Z(t) is the position operator in the Heisenberg representation, while 5(E — H) is the spec- 
tral measure of the SiNW Hamiltonian. Periodic boundary conditions are applied along the 
nanowire, the convergence being achieved for supercell lengths L ~ 500 nm. The real-space 
methodology of Ref. [19| has been adapted to this multi-orbital per site problem. In par- 



ticular, we have used the kernel polynomia 
from the Lanczos recursion coefficients 



method 



231 ] to compute the spectral quantities 



24J ]. This provides a more accurate description of 



the band edges than the usual continued fraction expansion [15|, |l9|, 125 ]. 
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FIG. 2: (color online). Mean free path of the holes (a) and electrons (b) as a function of energy, 
for L r = 0.54 nm (blue dashed line) and L r = 2.17 nm (red solid line). The electron mobility is 
also shown for L r = 2.17 nm (black dash-dotted line). 

RESULTS 

Electronic properties 

Let us first discuss the electronic structure of the SiNWs. The conduction band density 
of states (DoS) for the ideal and for two disordered SiNWs are shown in Fig. [T](b). In the 
ideal SiNW, the first two Van Hove singularities (VHs) arise from the [001] bulk conduction 



band minima [17[ and are split by the inter- valley couplings, while the others (above 1.7 
eV) arise from both [001] and {[100], [010]} minima. The DoS is markedly affected by the 
SRD. At L r = 0.54 nm, the lowest-lying VHs are shifted to higher energies, as a result of 
the increase of the average lateral confinement within the SiNW. However, with increasing 
L r , the conduction band edge moves to lower energy, while the DoS is steadily degraded, 
hardly showing any fine structure for L r > 2.17 nm. The lowest-lying electron states are 
indeed trapped deeper in energy in the largest sections of the nanowire [see Fig. (JT])a]. The 
extendedness of the electron wave functions will ultimately determine the transport regime 



Elastic mean free paths and charge mobilities 

The mean free path of the electrons and holes are plotted as a function of energy in Fig. [2, 
for the same two L r 's as in Fig. [TJ The main features of the underlying band structure still 
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FIG. 3: (color online). Maximum mobility for the electrons (red triangles) and holes (circles: first 
peak, squares: second peak) as a function of L r . Inset: Maximum mean free path. 

show up at L r = 0.54 nm. Indeed, the electron mean free path reaches its maximum 
(£ e ~ 70 nm) between the first two VHs (single subband transport), shows a dip at the edge 
of the second subband, then further decreases above E ~ 1.7 eV due to enhanced interband 
scattering. The hole mean free path also exhibits a very sharp peak (£ e ~ 110 nm) in the first 
subband, but becomes very short in the dense lower-lying subbands. However, at L r = 2.17 
nm, the fine structure of £ e can not be so easily related to the band structure of the ideal 
nanowire. The mean free path is almost reduced by half in the first electron "subband", 
while the peak on the valence band side is nearly five times smaller and is superseded by a 
broad feature at lower energies. 

The charge carrier mobility fi is another key quantity for assessing device performances. 
It can be related to the Kubo conductivity cr(E) = n(E)jjJE}e, where n(E) is the charge 



density (per unit of length) and e is the elementary charge |l9l.l20l|. Alternatively fi(E) can be 
written as er e {E)/m* , where m* is the effective mass of the carriers. These two expressions 
are fully equivalent under the assumptions of Drude's phenomenology with parabolic bands. 
However, the latter is impractical in SiNWs, because the bands have different effective 
masses (and might show significant non-parabolicity) , and also because the SRD induces 
large changes in the electronic structure as evidenced previously. /J*(E) is shown in Fig. [2]^b) 
for L r = 2.17 nm, and follows the same trends as £ e . 

The maximum mean free paths £™ ax and mobilities [i ma x are reported as a function of 
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L r in Fig. [3] for both kinds of charge carriers. The heights of the first, sharp peak and of 
the second, broad peak at lower energies are both given for the holes [arrows in Fig. 12(b)]. 
The mobility is always maximum at the first peak, although the mean free path can be 
larger within the second one. The values of n ma x, that are ranging from several hundreds to 
thousands of ran 2 V -1 s -1 , turn out to be in fair agreement with the experimental estimates for 
the most performant undoped semiconducting nanowires Q. The mobility can nonetheless 
be much smaller around the band edges, as evidenced in Fig. [2(b). Additionally, i™ ax and 
Umax show a minimum as a function of L r . This suggests that electrons and holes are less 
sensitive to short length scale fluctuations of the SRD profile (at least in the first conduction 
and valence subbands) [27]. The carriers are much more efficiently scattered by the SRD 
at intermediate values of L r . Finally, the surface of the nanowires becomes locally smooth 
again as L r is further increased, which enhances the mobility of the charge carriers far enough 
from the conduction and valence band edges. 



Conductance, quantum interference effects and localization regime 

For an in-depth analysis of transport properties in the coherent regime, the effects of 
quantum interferences effects on conductance scaling need to be analyzed in details. This 



can be conveniently achieved by using the Landauer-Biittiker transport framework [20(. This 
method is additionally well suited for describing transport between contacts connected to 
the nanowire, which becomes important when charge injection properties start to play an 
important. Properties such as transmission and local density of states are extracted from 
the Green's function and the self-energies of the contacts. 

r n (E)=i[Z n (E)-Zl(E)}, (4) 

T nm (E) = Tr[T n {E)GT m G\E)\, (5) 

where Y, n (E) is the self-energy of contact n and T nm (E) is the transmission between contact n 
and m. As the Green's function is calculated in segments using a Green's function decimation 



technique, 30f] it is an 0(R 6 ) method with regard to the nanowire radius, limiting practical 



calculations to R S, 1 nm. The contacts are modelled as semi-infinite nanowire leads with no 
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E= 1.475 eV 



E= 1.550 eV 
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FIG. 4: (color online), (a) Local density of states in a short nanowire with length L ~ 30 nm at 
energies i£ = 1.475 eV and E = 1.55 eV [arrows in (b)]. (b) Localization length £ as a function of 
energy (L r = 2.17 nm). (c) Linear fit (red solid line) to {\ng)(L) (blue crosses) at E = 1.55 eV. 



rou ghn ess and a radius of R = Rq + 0.2 nm. A standard decimation technique is employed 



15 



301 ] to compute the transmission through the SiNWs. 



In absence of any SDR, the transport regime remains ballistic and the resulting computed 
Landauer conductance G is quantized G ~ N±G |36|, with G = 2e 2 /h being the quantum 
of conductance, and N± gives the number of conducting channels. Conductance quanti- 
zation has been experimentally observed in narrow constrictions (quantum point contacts) 



made on silicon |3l| and short silicon wires |7J, |32| . However, most of the available results 
performed on longer and wider lithographically defined Silicon nanowires evidence the strong 
contribution of disorder effects (due to wire width fluctuations or charged impurities) that 
produce quantum interferences, localization [33j and charging effects [34j. 135]. 

In this work, the effect of electron-electron interactions is neglected, and the focus is put 
on the transition from weak to strong localization. Indeed, in the absence of inelastic scat- 
tering, quantum interferences build up beyond the diffusive regime, lead to the localization 



of all wave functions in the zero-temperature limit III, |28|, |29j. The localization length £(E) 
can be extracted from a scaling analysis of the conductance through (lng(E)) ~ —2L/£(E), 
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FIG. 5: (color online). £/(2£ e ) (red squares) and N± + 1 (green solid line) as a function of energy. 

where g = G/Gq and L is the length of the wire [Fig. life)]. £(E) is plotted as a function 
of the conduction band energy in Fig. Efb) (L r = 2.17 nm). The logarithm of g was aver- 
aged (for each L) over 150 random SRD profiles is provided in the inset. £ ranges from a 
few nanometers close to the gap up to ~ 500 nm at higher energies. The electrons indeed 
tend to localize in the largest sections of the nanowires near the conduction band edge [see 
Fig. ([T])a]. As an other illustration, the local density of states is shown at two conduction 
band energies in Fig. Ufa), for a short nanowire with length L ~ 30 nm. The corresponding 
localization lengths are £(1.475 eV) = 13.5 nm and £(1.55 eV) = 114.5 nm. 

In weakly disordered quasi-l D sy stems, the fundamental length scales £ and l e are shown 



to scale as £ ~ 2(iVj_ + l)£ e 28, 37]. This relation was first established by Thouless [28] for 
strictly ID systems, and further generalize with random matrix theory to weakly disordered 



29 
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Recently, it was studied 
In our situation, the 



quasi-lD systems with larger number of conducting channels 
numerically for chemically doped disordered carbon nanotubes 
identification of a well-defined channel structure is however hindered in the SiN Ws by the 
SRD-induced changes of the electronic structure. The comparison between £/(2£ e ) and 
N±_ + 1 (deduced from the band structure of the ideal nanowire) is nonetheless instructive 
[Fig. [5]. As evidenced by our calculations, £/(2£ e ) roughly scales as N± + 1 close to the 
conduction band edge, but the situation becomes much more complex at higher energies 
at which strong difference occur. Similar discrepancies were recently reported in simplified 



models of disordered quantum wires, but in the presence of a magnetic field 12]. 
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CONCLUSION 



In conclusion, some key transport length scales in disordered semiconducting nanowires 
have been investigated using an optimized real space order N method combined with a recur- 
sive Green's function-based Landauer-Biittiker approach. The effectiveness of an atomistic- 
based surface roughness profile in limiting ballistic transport has been demonstrated, and 
the trends in the energy-dependent electron and hole mobilities, mean free paths and local- 
ization lengths have been discussed at a quantitative level. The limitations of the Thouless 
relationship in such complex disordered systems have also been pointed out. Studies focusing 



on t 



17 



re role of nanowires orientation, diameter, and other kind of disorder such as dopants 



391 ]. surface defects, or traps in the oxide, deserve further consideration. Additionally, 



beyond the intrinsic effects of these various sources of scattering on transport length scales, 
their impact in SiNWs-based field effect devices should be investigated [40]. 
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